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We consider general mixtures of isocurvature and adiabatic cosmological perturbations. With a 
minimal assumption set consisting of the linearized Einstein equations and a primordial perfect fluid 
we derive the second-order action and its curvature variables. We also allow for varying equation 
of state and speed of sound profiles. The derivation is therefore carried out at the same level of 
generality that has been achieved for adiabatic modes before. As a result we find a new conserved 
super-horizon quantity and relate it to the adiabatically conserved curvature perturbation. Finally 
we demonstrate how the formalism can be applied by considering a Chaplygin gas-like primordial 
matter model, finding two scale-invariant solutions for structure formation. 
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I. INTRODUCTION 

Cosmology, as every other field in science, comes 
equipped with its set of paradigms. Here we wish to 
break with the assumption that cosmological perturba- 
tions are adiabatic. Departure from this claim has of 
course already been discussed and considered in many 
forms in the literature (see e.g. PHI])- The novelty of 
this paper then is to develop the cosmological perturba- 
tion formalism in a relativistic context for any mixture 
of entropy and adiabatic perturbations in the same gen- 
erality that has been achieved for purely adiabatic per- 
turbations before (see e.g. Q). As such we will carry 
out our derivations employing a minimal assumption set 
only consisting of 

• the linearized Einstein equations and 

• modelling primordial matter as a perfect fluid. 

Importantly we therefore do not constrain the equation 
of state or speed of sound profile of the primordial fluid, 
which are consequently allowed to vary freely. 

In a nutshell our philosophy is therefore to say as much 
as possible about non-adiabatic perturbations whilst 
specifying as little as possible w.r.t. the nature of the 
underlying primordial matter. The resulting formalism 
will allow straightforward computation and comparison 
of e.g. power spectra produced by any model satisfying 
our assumption set. In contrast, previous literature has 
focussed on considering isocurvature perturbations pro- 
duced by more tightly constrained matter in great detail 
(e.g. matter-radiation and other fluid mixtures or 
multiple scalar fields Q)- 

Observationally, purely adiabatic modes are still com- 
patible with current experimental data, but near-future 
experiments such as Planck should soon tighten con- 
straints, possibly establishing the presence of entropic 
modes. In fact there might already be positive evidence 
suggesting their presence (see e.g. 

The paper is organized as follows. In section|H]we give 
a precise definition of what we mean by "non-adiabatic" 
perturbations and give explicit examples to show how 



isocurvature modes can arise from the interaction of com- 
ponents in an overall perfect fluid. In Mil we then extract 
equations of motion from the linearized Einstein equa- 
tions for a perfect fluid. These are developed into a 
second-order action for fluctuations and the associated 
curvature variables in the adiabatic case in IIVI Here we 
also show that single scalar fields generically have van- 
ishing non-adiabatic perturbations on large scales. In 
section [V] we then present the main results w.r.t. the 
formalism developed in this paper: A second-order ac- 
tion for fluctuations and associated curvature variables 
in the non-adiabatic case. We pay special attention to 
the asymptotic limits of the solution, finding a new con- 
served super-horizon charge. 

Having mentioned some observational constraints im- 
posed on primordial matter models in IVI1 in section IVIII 
we provide an application of the formalism developed. 
More specifically we investigate a Chaplygin gas-like pri- 
mordial matter model, finding two solutions capable of 
producing scale-invariant power spectra and resolving the 
horizon problem. Finally we look at the issue of growing 
and decaying modes in I Villi before concluding in section 

El 



II. (NON-)ADIABATICITY 

A. The hydrodynamical setup 

Let us begin by spelling out what exactly we mean by 
"non-adiabatic perturbations" . In staying with our phi- 
losophy of leaving the exact nature of the primordial mat- 
ter in question unspecified, the formalism developed here 
will be at the level of hydrodynamics and we will thus 
characterize the system thermodynamically by its pres- 
sure (p), energy density (e) and entropy (S) variables. 
When Taylor-expanding perturbations in the pressure in 
terms of the other variables we find 
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where we have denoted the entropic perturbation term 
as 8p en t and have defined the adiabatic speed of sound 
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where ' is a derivative w.r.t. proper time t. 

Our demarcation criterion in distinguishing non- 
adiabatic from adiabatic perturbations will therefore sim- 
ply be Sp en t 0. We therefore group together purely 
isocurvature perturbations (for which Se = and hence 
Sp = Sp ent ) and arbitrary mixtures of isocurvature and 
adiabatic modes, using " non-adiabatic " as an umbrella 
term to describe all overall pressure perturbations with 
non-zero Sp en t- It is worth emphasizing that this means 
we take "non-adiabatic" to literally mean "not adia- 
batic", and not only as a reference to purely isocurvature 
modes. 

For notational convenience we also define an "effective" 
speed of sound c 2 s here, which relates pressure and energy 
density perturbations and reduces to c s 2 when evaluated 
at constant entropy, i.e. in the adiabatic limit 
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B. A perfect fluid 



(4) 



The primary assumption going into our derivation will 
be treating the overall primordial matter as an effective 
perfect fluid (as it arises naturally for FRW universes). 
Consequently our primordial matter stress-energy tensor 
will be of the following form 







-p)u a up -p5 a p. 



(5) 



u here refers to the potential 4-velocity. 

Specific single fluid/field examples of matter with such 
a perfect fluid T a @ include hydrodynamical fluids as well 
as scalar fields cj> with Lagrangian = y/—gp(X, </>), 
where the kinetic term is X = \g a $ § _p . As such this 
incorporates canonical scalar fields with potential V((f>) 
as well as those with non-canonical kinetic terms (e.g. 
as encountered in k-inflation [SI). For a canonical scalar 
field (j) we obtain (JSJ with (see Q) 



P = x - v{4>) 



e = X + V(4>). 



(6) 



More generally, however, we can view such an overall 
effective T a p as representing the sum of stress-energy 
tensors corresponding to (possibly multiple) components 
and their interactions. It is these interactions that can in 
principle give rise to non-adiabaticity in an overall perfect 
fluid. 

To illustrate the mechanism by which this can happen, 
consider a perfect fluid with the following two-component 
decomposition 

T a p = (e + p)u a up - P 8 a fj = Tx a p + T 2 V (7) 



then show that we can derive the following relation- 
ship independently of the nature of interactions between 
the two components 



P-\D = ^[p,~^D,\^[p 2 -^D, 
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Here P and D are measures of the overall fractional pres- 
sure and energy density perturbations, where subindices 
i denote corresponding quantities for individual compo- 
nents i. 
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From our discussion above it is clear that = c?„ 



(9) 



(10) 



and 



hence that — c s 2 = % in the adiabatic case. The first 

D S £ 

two terms in therefore vanish if the individual compo- 
nents behave adiabatically. However, the third term is a 
measure of the interaction between the two components 
and is non-zero in general. From this we can immediately 
conclude that for the overall fluid = c 2 es ^ | = c s 2 . 
In other words we can obtain non-adiabatic fluctuations 
in the overall perfect fluid in question. A closer look at 
the [] bracket term in also shows that departure from 
adiabaticity is directly linked to having distinct speeds 
of sound for individual components. 

As an explicit example we here mention the familiar 
case of the Chaplygin gas, whose popularity stems from 
potentially providing an effective single fluid description 
of dark energy and dark matter. A Chaplygin gas has 
equation of state p = ^. In its simplest forms A is a 
constant, however here we will allow fluctuations in A as 
well, following once again. One then finds that the 
perfect fluid Chaplygin gas can be viewed as composed 
of another Chaplygin gas interacting with a pressureless 
fluid. Overall we obtain 



P 



,D 



t 1 = 
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SA 
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(11) 



(12) 



For a non-zero p, we therefore again find c 2 s ^ c s 2 and 
have produced a perfect fluid model with non-adiabatic 
perturbations. 



III. PERTURBATION EQUATIONS 

Here we derive the the equations of motion associated 
with relativistic primordial perturbations from the lin- 
earized Einstein equations. Establishing our notation we 
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write the Friedmann equations as 

8ttG 



U 2 = 



a 2 e 



H 2 -H =4irGa 2 (£+ P ) = G, 



(13) 



(14) 



where e and p represent the overall background energy 
density and pressure respectively. We also have the con- 
formal Hubble factor % = a' /a, a(i]) is the scale factor 
of the FRW metric and Q has been defined in order to 
keep notation concise in what follows. Here ' denotes 
differentiation w.r.t conformal time 77, which is defined 



drj = 



dt 
a 



(15) 



The linearized perturbed Einstein equations for a per- 
fect fluid T a p © are 



A$-3"H($ + H$) = 4irGa 2 5e, 



(16) 



In the adiabatic case the equations of motion therefore 
simplify considerably. However, since we will eventually 
want to drop the requirement of adiabaticity, we proceed 

differently here. We combine the Su^ and 6u\\ terms by 

introducing the new variable £ = duuf, requiring 
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Importantly / is therefore only defined up to its frac- 
tional variation here. Comparison with (|23p shows that 
/ provides a generic measure of the departure from adia- 
baticity. If we furthermore define £ = a<3>, we have com- 
pactified our equations of motion (|2ip and (|2"2"j) to the 
following form 



? r 



c = /(i + 



(25) 



(26) 



(0$) i = 4irGa 2 (e +p)5un = QSi 



(17) 



+ 3"H$' + (2H + H 2 )$ = 4irGa 2 5p. (18) 

Here A is the Laplacian, $ is the Newtonian potential 
and the stress-energy fluctuations are evaluated in the 
longitudinal gauge (denoted by an overbar) ; we refer the 
reader to [1,0 for notation and further explanations. We 
use equation (fT7|) to solve for $ and <f> : 



Q — 

$ = -bu\\ -U<f> 
a 11 



(19) 



$ = — bun H — (- 



p 



- 2H)6u\\ + G<$> (20) 



a " a e + p 
Combining equations (fT5|) and (|T5)) we get 

d>" + 3(1 + cl s )U$ + (2H + (1 + 3c 2 es )H 2 - c 2 s A)$ = 

,(21) 

Using P^|) and ([2T7| to substitute in for $' and we 
also derive the following equation for 8u\\ 



5v,n + (5zt|| ( 



(e+ P ) 



3c 2 es H) = (1 + 



-)(o$) (22) 



Equations f|2 1 [) and (|22|) are effectively the equations of 
motion for the system under consideration, governing the 
Newtonian potential $ and the matter 4-velocity u re- 
spectively. They are also the precursors to the so-called 
u- and D-equations we encounter in the literature. 

Noticeably the Sun term in (|2"2"]l vanishes for adiabatic 
perturbations, since 



P 



Zc 2 U = 0. 



(23) 



Using the methods outlined in Appendix A we can now 
also write down an action for this system 



'v' 2 



ci s A)v 2 + —v 2 )dr)d 6 x, 



where v = zC, and we have defined 



z 2 = 



go 



f(Q + 4AY 



(27) 



(28) 



where 0(A) is some time- independent operator (see ap- 
pendix A for details). Whilst (|27[) is a perfectly valid 
action, however, it is representationally rather opaque. 
For instance, multiple terms here depend on A and hence 
(after Fourier-transformation) on the wave mode k. This 
obscures physically significant features such as the behav- 
ior of long- and short-wavelength modes, so that we will 
find it helpful to recast (|27l) into a form where such fea- 
tures as well as the relative contributions of entropic and 
adiabatic modes become more apparent. The next two 
sections are dedicated to this task, first in the adiabatic 
limit and then in the non-adiabatic case. 



IV. THE ADIABATIC CASE 
A. General solution 

Here we are considering the limit where Sp en t = and 
hence /' = 0. Using the techniques outlined in Appendix 
A we transform the equations of motion ((25)) and f|26[) by 
switching to new variables 



Cad — 



Cad 



a a 



(29) 
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In terms of these new variables we obtain equations of 
motion 

n 2 



Ga 2 

The associated action is 



Cad — -TJjCad Cad ~ T r ^l C ^ '^Cad- (30) 

ti l Ga 



(v' 2 + c s 2 Av 2 + — v 2 )drid 3 x, (31) 

z 

where v is the familiar canonical quantization variable 

v = z( ad and 



z 2 oc 



Qa 2 Q 
n 2 c s 2 A' 



(32) 



Identifying the time-independent operator O with the 
Laplacian, O = A, we have reproduced the results ob- 
tained in the literature for adiabatic fluctuations Q. 
z 2 is independent of k and the representational advan- 
tages discussed above are restored. Furthermore C is the 
well-known curvature perturbation, which freezes out in 
the long- wavelength limit, where we can ignore the c 2 
term. This feature is easily visible, once we've varied S 
w.r.t v and Fourier-transformed the resulting so-called 
-^-equation, obtaining 



2;„2 



v + (eg k 



0. 



(33) 



B. A single scalar field 



Perhaps the most popular primordial matter candidate 
in the literature is a single scalar field. We will here sum- 
marize the argument given in (Tlj . showing that a single 
p(X, 4>) scalar field only produces adiabatic perturbations 
on large scales, as identically Sp ent = for any such field 
then. For we can write the entropic pressure Sp en t as 



Sp e 



P,<t> ( l ' 
P,x (l 



2 ) - 2c 2 Xp. x 
2 ) - 2c 2 Xp 



XX 



8X. 



(34) 



In addition we can derive the following constraint equa- 
tion on large (super-horizon) scales 



\l>..\ - 1 r I'. A v I SX 

[4> 2 P,x<t>- p,<j, + 3Hp 



',x(, 



0. 



(35) 



Substituting this into the background field equations 
leads to SX = <j>S4>. Finally, substituting this and the 
following expression for c 2 



P,X(t)+P,4> 



(36) 



P,x4> - p,4> + P,xx4> 2 § + P,x^ 2 
into ([M]) . we find that 

Spent = 0. 

A single p(x, 4>) scalar field therefore only produces adi- 
abatic perturbations on large scales, i.e. c 2 s = c 2 . 



(37) 



V. THE NON- ADIABATIC CASE 

A. General solution 

What happens if we allow Sp en t to depart from 0? Let 
us firstly recast the equations of motion into a form, 
where we separate out dependence on adiabatic and 
isocurvature modes. In terms of the adiabatic variables 
C, a d and ^ a d we now obtain 



Go 2 
H 2 



£ad ~ 1/2 Cad 



n 



Cad ~ n y C es^Cad + . ( ? Cod Ca 



(38) 



(39) 



This explicitly reproduces the adiabatic equations in the 
/' = limit and we can identify the last two terms in 
(|39f as "non-adiabatic" correction terms. Their depen- 
dence on / // shows that they vanish identically in the 
adiabatic limit. We have therefore separated out adia- 
batic and isocurvature modes at the level of the system's 
equations of motion. 

In order to extract information about the behavior of 
asymptotic wave modes, it is, however, helpful to com- 
bine those separate contributions into a more compact 
form. Following this agenda, equations (|38f and ((39)) 
can be further compactified by iteratively applying sub- 
stitution schemes in analogy to the derivation of equa- 
tions ([2"5"]l and (|2l)]). After a series of such steps we have 
mapped variables Cad and Cad to 



= Ai/Cad - ^Cad 



p ; 



(40) 



Here the new substitution functions v and [i are defined 
via coupled differential equations |39j | 



Hf 



Ga 2 

n 2 f l 



(41) 



With this new set of functions at our disposal we can 
finally write down the resulting equations of motion for 
the system 



P = 



Ga 2 1 
H 2 f(i 2 



n 2 f 2 2 



Ga 2 



P c es Ap, 



(42) 



(43) 



Using techniques from Appendix A we can turn these 
equations of motion into an action 



S = I {v 2 + c 2 ;s Av 2 + —v 2 )dr)d 3 x, (44) 
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where the associated curvature variables v and z are given 
by 



z 2 oc 



Go 2 



(45) 



Here we have once again identified O with the Laplacian 
A. 

With these variables the equation of motion for v is 
formally identical to that derived for the adiabatic case 



v — c 2 Av v = 0. 



z 
z 



(46) 



All the effects of non-adiabaticity have been absorbed 
into c 2 s and new expressions for z 2 and v (|45[) (and con- 
sequently a new time-dependent mass term z jz). 



B. Asymptotic limits 

Let us first consider the short- wavelength limit. Here 
the asymptotic solution to the fluctuation equations re- 
mains unaltered. To see this we Fourier-transform the v 
equation, obtaining 



(ctk 2 



0. 



(47) 



The two v-terms directly yield the short- and long- 
wavelength plane wave perturbations respectively. In 
the limit where the pressure term dominates (and hence 
c 2 k 2 » \z /z\) we can use the WKB approximation to 
give the short-wavelength solution 



1 



\J c es k 



exp(±ik / c es drf) 



(48) 



The functional form of the asymptotic short-wavelength 
limit is therefore independent of z / z and hence remains 
unaffected by the introduction of non-adiabatic modes. 

The long-wavelength solution is obtained by consid- 
ering the opposite limit when the "time-variable mass" 
term dominates the pressure term (i.e. when c 2 k 2 << 
\z /z\). The introduction of entropy modes therefore 
plays a significant role. 



v = C 1 {k)z + C 2 {k)z 



drj 



+ 0((k v 2 )), (49) 



where C±(k) and C2{k) are time-independent functions 
of the wavenumber k. Typically the first term in this 
expression represents the growing and therefore domi- 
nant mode allowing us to ignore the second term. We 
will discuss the effects of a dominant C 2 term in IVIIII 
For a dominant growing C\ mode, however, one imme- 
diately obtains a conserved quantity on super-horizon 
scales, i.e. in the long-wavelength limit. Following our 
notation above this is given by 6, for v — zO (see (145(1 ^1 . 



thus resulting in a simple conservation equation 9 = 0. 
This reduces to the familiar C, ad — in the adiabatic 
limit. The "frozen-in" quantity outside the horizon, still 
formally given by v/z, therefore generalizes in the way 
presented. 

Relating the new expression for 6 to the adiabatically 
conserved Cad one finds 



Cad 



(v£ad) 



(— + -j)Cad- 

M J 



(50) 



The right-hand side of this equation can now be in- 
terpreted as an entropy producing source-term caus- 
ing a deviation from the previously conserved Cad- 
This has important physical consequences. For exam- 
ple non-linearities could now continue to grow outside 
the horizon, potentially generating large levels of non- 
Gaussianity (l2j ] . 

Finally we can also obtain a general expression for hori- 
zon crossing modes, when considering the point where the 
short- and long- wavelength solutions "meet". Modes 
crossing from inside the Hubble radius to the gravity 
dominated region outside are then subject to the require- 
ment that 



(51) 



(*•) = "FT > 0- 



If this inequality is not satisfied, then the modes consid- 
ered are "moving" in the opposite direction, i.e. crossing 
the Hubble radius from the outside in. 



VI. COSMOLOGICAL CONSTRAINTS 

The formalism developed in the previous sections is 
completely general, as long as our minimal set of assump- 
tions is met. In particular we have not constrained the 
equation of state of the underlying fluid or the associated 
behavior of the scale factor. Inflating or non-inflating 
(ii > or < 0) as well as expanding and contracting (a > or 
< 0) solutions are still all on the table. Here we will com- 
pile a simple checklist to classify solutions arising from 
any particular model we might consider. 

• Is the model capable of resolving the horizon prob- 
lem? 

• Is an "expanding" solution with a > possible? 

• Can a scale-invariant power spectrum of pertur- 
bations be produced? The observed near scale- 
invariance of fluctuations [l3| therefore completes 
our check-list with an observational criterion. 

Let us briefly expand on the items in this list. The hori- 
zon problem is posed by the fact that the scales we now 
observe are initially causally disconnected according to 
the unreformed Big Bang model yet are also stunningly 
uniform. In other words, the remarkably homogeneous 
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and isotropic energy density distribution we observe in- 
side the present horizon scale rfo ~ 10 28 cm appears to 
require a strongly fine-tuned initial matter distribution 
in the absence of a physical mechanism generating such 
uniformity. Schematically a "resolution of the horizon 
problem" involves a temporary reversal of the kinematics 
realized in the unreformed Big Bang model. Specifically 
this requires the realization of a phase in the early uni- 
verse when modes cross from the region inside the Hub- 
ble radius (where they are ruled by causal micro-physics) 
to the region outside (where they become dominated by 
gravity). With respect to the formalism this means that 
((STj) needs to be satisfied by fc-modes during this phase. 

In this context inflation is often invoked where 
initial vacuum quantum fluctuations inevitably present 
while the modes are inside the Hubble radius are 
stretched to super-horizon scales during a phase with ac- 
celerating scale factor d(t) > 0. Typically this is achieved 
via the introduction of one or multiple scalar fields. Ex- 
amples include slow-roll inflation, where restrictions are 
placed on the form of the scalar field's potential V{<j>) (see 
e.g. UM) 

and a variety of models with non-canonical 
kinetic terms (see e.g. [t|[l7j and references therein for k- 
inflation and DBI-inflation respectively). Many alterna- 
tives exist, however: Cyclic or ekpyrotic scenarios [l~§.ll~9l] 
propose a contracting pre-Big Bang phase " causally con- 
necting" fluctuation modes. String gas cosmology [2(ij 
invokes thermal fluctuations that exit the Hubble radius 
during an early quasi-static Hagedo rn p hase. Varying 
speed of light (or sound) frameworks |2lTl23l | spread per- 
turbations via a much larger speed of sound in the early 
universe, thus allowing all observed modes to originate 
from causally connected regions. Combinations of all of 
the above (see for example [12|]) have also been proposed. 
More recently yet another potential resolution has been 
pointed out within Hof ava-Lifshitz gravity [24l428l | . Here 
the effective speed of light (and hence also the speed of 
sound) diverges in the ultraviolet regime, thus resolving 
the horizon problem in a way conceptually analogous to 
varying speed of light models. 

The majority of these proposed "resolutions" to the 
horizon problem are realized within an expanding phase 
of the universe (a > 0). However, there are exceptions; 
e.g. ekpyrotic scenarios or some of the models consid- 
ered by us in [29j]. Here the evolution is reversed as a 
contracting universe is getting denser and hotter in time 
rather than diluting and cooling. Features such as scale- 
invariant spectra are set up within a contracting phase 
in these models and the horizon problem is solved by fiat 
as causal connectivity is invariably produced during the 
contraction. In section I Villi we point out some of the 
issues one must address when employing a contracting 
phase to generate fluctuations. 

One should note that there also remains a gauge prob- 
lem in dealing with perturbations here. This is due to the 
fact that varying w models give rise to different scales 
for metric and matter variables. We will not replicate 
the discussion of this topic here, but refer the interested 



reader to [Hj]. 

VII. A MODEL EXAMPLE: THE CHAPLYGIN 
GAS REVISITED 

We will now provide an example of how the for- 
malism developed can be used to extract features for 
non-adiabatic perturbations around some specified back- 
ground. The particular example we choose is an exten- 
sion of the model discussed in (2||. Here the equation 
of state w and the effective speed of sound c 2 ^ exhibit a 
power-law dependence on the energy density s. In prin- 
ciple w and c 2 s of course do not have to be functions of 
e only, but for illustrative purposes we will restrict our- 
selves to this simple model here. Other models exhibiting 
such a dependence include e.g. the Chaplygin gas and 
its modifications [30j ] . but also intermediate inflationary 
models [3l|, 132J . In fact the latter leads to a second so- 
lution for scale-invariance in the inflationary setting (the 
other solution being slow-roll inflation (l6l. l33j). 

A. The background solution 

We consider models with equation of state w = p/e 
and effective speed of sound c 2 s given by 

w + Wqcx e 2fi <4 cx e 2a (52) 

In the adiabatic limit c 2 s = c 2 , which corresponds to 
a — (3 here. Specifying the equation of state in the adi- 
abatic case will therefore automatically fix c 2 as well, 
so that we can then completely describe the solution via 
one parameter. In contrast, in the non-adiabatic case c 2 s 
in principle becomes a free variable, so that we need at 
least two parameters to specify our model. 

We will discuss solutions with positive and non-zero (3 
and positive a here. For other choices of a, /3 we refer to 
appendix B. At high density these models all display the 
same (power-law) behavior. A high 10 phase then exits 
into a constant w, low density phase which need not be 
inflation. The corresponding low energy equation of state 
is given by wq- We can also consider contracting models 
with these equations of state: the obvious generalization 
of cyclic models, where a constant high w is invoked. 

Solving for e we get 

eoct^+h, (53) 

where we have ignored the low-energy contribution Wq. 
The energy density therefore changes like a power-law in 
t and diverges as t — > 0. As such we have 

-40 „ -4a 

w cx <f es cx t 1 ^ (54) 

Furthermore we get the following expression for the scale 
factor a(t) 

a(t) = a e K « tT ^ « a . (55) 
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where the approximate equality holds in the regime we 
shall be interested in, namely for early times t « 1. The 
universe therefore appears to be loitering as a function 
of time (a ~ ao) and a does not vanish at t = 0. This 
means we can effectively use conformal and proper time 
interchangeably, as they only differ by a constant factor: 
i] = t/ao when t « 1. Nevertheless the fact that the 
universe's energy density diverges as t — > shows that 
we do still have a Big Bang singularity. Completing our 
background solution we find that the conformal Hubble 
factor behaves as 

H oc at 1 ^ k aot 1 ^ . (56) 

B. Perturbations and Power Spectra 

Having specified an equation of state w and a partic- 
ular profile for c 2 s we are now in a position to explicitly 
solve the fluctuation equation (|47D . In order to work out 
the associated power spectra of fluctuation modes we in 
principle need to know the full Bcssel function solution 
for v. However, we can in fact extract all the neces- 
sary information from the asymptotic solutions to the v 
equation ([35)1 and flU?)) , When the growing "Ci"-mode 
dominates, we can simplify these expressions for long- 
and short-wavelength modes respectively to the follow- 
ing form 

v oc Ci(k)z v cx — ^=r, (57) 

V c esk 

where we have ignored a phase for the short wavelength 
limit. We can now determine Ci(fc) by requiring these 
two asymptotic solutions to match when the z;-term van- 
ishes, i.e. at the point of horizon-crossing. Consequently 
we "glue" solutions at the point where 

cl s k 2 = Z -. (58) 

The scalar power spectrum for v, \Ve(k)\ 2 , is given by 

\V 9 (k)\ 2 ock 3 \9\ 2 , (59) 
which we can express at the point of gluing as 

|p e (fe)| 2 ocfc 3 |C 1 (fe)| 2 c<4^ = ^^- (60) 

For a scale-invariant solution \Pe{k)\ 2 has to be dimen- 
sionless. In other words |Ci(£;)| 2 oc k~ 3 . Furthermore the 
dimensionless amplitude of the power spectrum is experi- 
mentally well constrained. We should also point out that 
(|60[) is valid for modes exiting the horizon at a primordial 
stage regardless of the form of the long- wavelength solu- 
tion for v. However, for it to still be a valid expression for 
modes re-entering the horizon at a later stage we require 
a dominant C\ mode (see section [Villi for details). 



Substituting from the background solution into the 
definition of z 2 ([15)1 and differentiating, we also find that 

/ 4(q-/3) j /' ti a 

— oc — —t -r — 2 h2 — . (61) 

z 1 + 4/3 / fx a K ' 

This is interesting, as it contains information about which 
terms will dominate the evolution. Note that from ([56)) 
we can see that the first term will generically suppress 
the a I a = % contribution, unless a = /3, i.e. the evo- 
lution is adiabatic. As such we can typically ignore the 
contribution of the T-L term in non-adiabatic cases. This 
is intriguing, given that it is the % term that provides 
the dominant contribution in familiar (adiabatic) infla- 
tionary models where z 2 oc a 2 . With an approximately 
constant ao, however, we still obtain z jz oc r\~ 2 when 
only the t~ l term is significant, which is reminiscent of 
the analogous z / z oc ri~ 2 expression we frequently find 
in inflationary models [f 21 ]. 

C. Adiabatic scale invariance 

As can be seen from ([52| . adiabaticity corresponds to 
a = P (with appropriately tuned proportionality con- 
stants). We here recover the familiar z 2 oc a 2 , which re- 
sults in the following expression for fc-modes in the early 
time limit 

(k 2 )' oc (H 2 )' < 0. (62) 

No expanding resolution of the horizon problem is there- 
fore possible here, as modes are crossing over from out- 
side the Hubble radius to the inside region. This is the 
opposite behavior to the one necessary to realize the de- 
sired reversal of the unreformed Big Bang model kine- 
matics. In order to address the horizon problem this 
model consequently has to be implemented within a con- 
tracting phase, as e.g. invoked in cyclic scenarios. For 
in such a phase the evolution is reversed and incoming 
modes are mapped into outgoing ones and vice versa. 

Applying the gluing techniques discussed above we can 
now compute the scalar power spectrum for adiabatic 
perturbations at the point of horizon-crossing. We find 

\V e (k)\ 2 ocfw, (63) 

straightforwardly giving us a criterion for scale invariance 
by requiring the spectrum to be k- and hence, at horizon- 
crossing, time-independent 

P = a = l. (64) 

This result is in agreement with the one obtained in ^20j j . 
We therefore recover a scale- invariant solution which fur- 
thermore is capable of solving the horizon problem in a 
collapsing phase. We should note that by simultaneously 
requiring adiabaticity and (3 > we have only explored 
time-varying c s 2 solutions here. 
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D. Non-adiabatic scale invariance I 

Let us first consider non-adiabaticity with a > (3. 
Comparing with (I61[) we then find 

— oc oc i^FT 1 (65) 
z J 

As such / will exhibit an exponential dependence. We 
again emphasize that the usually important contributions 
from the scale factor are negligible here. We are therefore 
in a radically different regime from the standard infla- 
tionary z I z oc a I a scenario. With this the fluctuations 
are fully specified and, after computing the time- variable 

mass term — to first order, we find 

2 )' cx (TW^)' < 0. (66) 

Analogously to our reasoning for the adiabatic case, we 
can thus conclude that no expanding resolution of the 
horizon problem is possible here as well. In fact this con- 
clusion is not an artefact of choosing a power-law solution 
for c 2 s , but will obtain for a much wider class of c 2 s pro- 
files, e.g. if an exponential trial solution c 2 s oc e kt is 
chosen for some constants c and k. In the class of scenar- 
ios considered here it would therefore require a change in 
the equation of state to overcome the solution's "affinity" 
towards collapsing models. For, whilst an initially large 
and subsequently decaying c 2 s has the tendency to re- 
verse the horizon crossing behavior of k-modes, the time 
variable mass term z j z counteracts this here. 

Moving on to computing the scalar power spectrum we 
obtain 

„ -2q-2 
V f 1 + 4/3 

|7M£0| 2 (x^^cx^^. (67) 

z c es z 

However, the exponential dependence of the /-measure 
(|55|) here means that z 2 diverges as t — > 0. And hence 
1/z 2 vanishes as t — > 0, growing with time afterwards. 
From expression (|60l) we see that the amplitude of the 
power spectrum \Vg(k)\ 2 will consequently grow with 
time, also showing divergent behavior. Comparing this 
with (|66D we find that the result is an exponentially red- 
tilted power spectrum. Baring extreme fine-tuning, no 
(near)scale-invariant solution can therefore be produced 
here. 



E. Non-adiabatic scale invariance II 

We will now consider the non-adiabatic regime with 
P > a. Prima facie one would not expect this regime 
to succeed in stretching primordial fluctuations to super- 
horizon scales in an expanding phase. For whilst we can 
still have an initially large speed of sound here, its effects 
are offset by a more rapidly varying equation of state w. 
Nevertheless it might present an interesting possibility 



for contracting scenarios. Let us first back up the claim 
that no expanding resolution of the horizon problem is 
possible in this context though. 

Regarding the formalism we find that the "non- 
adiabatic" functions / and fi are given by 

focQ fiocQ- 1 (68) 

For j3 > a > we consequently get horizon-crossing 
modes with 

(k 2 )' oc (f 2+ ™)' < 0. (69) 

This does indeed show that no expanding solution can 
solve the horizon problem in the proposed fashion here. 
One should point out that the reasoning for the case 
where a = fi, but adiabaticity fails due to non-matching 
proportionality constants, is entirely analogous to the 
considerations presented in this section and the same 
equations are reproduced. 

Is the j3 > a solution a candidate for a contracting 
solution to the horizon problem? Turning our attention 
to the scalar power spectrum we find that \Pe{k)\ 2 at 
horizon-crossing is given by 

|^(fc)| 2 cxtw, (70) 

For a scale-invariant solution we require ([70]) to not be 
t-depcndent, so that the condition for scale invariance is 
simply 

a = 0. (71) 

This corresponds to having a constant effective speed of 
sound c 2 s . We have therefore found a non-adiabatic, con- 
tracting and scale-invariant solution with constant c 2 s 
but arbitrary w oc e 2 ^ here. Interestingly it is straight- 
forward to fabricate blue-tilted spectra in the (3 > a > 
regime. This can be achieved by choosing positive, non- 
zero a as can be seen from ((TO)) . A red-tilted spectrum 
is not possible here, however. 

In summary we have found a family of equations of 
state w that satisfy our background solution with /3 > a, 
which can resolve the horizon problem and produce flat 
power spectra in a contracting phase. 

VIII. GROWING VS. DECAYING MODES 

We remind ourselves of the conservation equation for 
long-wavelength modes (pHJj) 

v = Ci(k)z + C 2 (k)z / J + OdkTj 2 )). (72) 

This will only give rise to a conserved quantity 9 if the 
C\ term is the dominant/growing mode. If the C2 mode 
is dominant, on the other hand, 8 is not conserved. A 
power spectrum scale-invariant in 9 upon horizon exit 
consequently does not guarantee that modes re-entering 
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the horizon at a later stage will be scale-invariant in 9 
anymore. Setting up a primordial scale-invariant power 
spectrum then no longer automatically justifies observ- 
ing a scale invariant spectrum today. Meeting the obser- 
vational constraints therefore requires a more elaborate 
setup than the simple gluing considerations we presented 
for a growing C% mode. As an example of such a scenario 
with dominant C 2 term we here refer to a single scalar 
field in a contracting phase ([H, HE|). Such setups fre- 
quently suffer from the problem of obtaining generically 
blue Cad fluctuations (for an exception see |19j). 

The solutions worked out in the previous sections must 
also be implemented in a collapsing phase if they are to 
resolve the horizon problem. It is thus natural to ask, 
whether we are also faced with a dominant C 2 mode 
here. For in constructing such a collapsing phase one 
effectively considers the mirror image of the post Big- 
Bang picture we are familiar with. As temporal behavior 
is reversed upon mirroring the dynamics around a sin- 
gularity at t = 0, growing modes can be mapped into 
decaying ones and vice versa. When implementing a so- 
lution in a contracting phase, one therefore needs to be 
careful to not ignore previously suppressed solutions that 
can become dominant in a collapsing phase. 

For the example model presented here the following 
picture emerges: In the adiabatic limit the C2 term is 
proportional to J -J oc \r)\. It consequently vanishes as 
t — > and only the C\ term stays relevant. 

When a > f3 > the C2 contribution exponentially 
tends to zero as t — > 0. This can be seen from the ex- 
ponential dependence of z 2 ([55]) . It therefore becomes 
suppressed and hence irrelevant in this limit. Thus we 
maintain a conserved super-horizon quantity 6 even in 
the case of a collapsing phase. 

When /3 > a > we find that the C\ and C 2 modes 
are proportional to 

z\ cx t 1 + 4 ' 3 z 2 I oc t 1 + 4 ' 3 . (73) 

J z 2 z 

From these expressions it again follows that the C\ mode 
does indeed dominate over the C2 mode and we therefore 
maintain a frozen-in 9 on super-horizon scales. 

We have therefore justified that the "gluing" procedure 
in our derivation does give us the correct power spectra 
both at horizon exit and re-entry for all cases consid- 
ered. In other words, we are protected from additional 
super-horizon fluctuations sourced by C2 modes. We will 
not discuss the issue of matching modes across the singu- 
larity here, but have worked under the assumption that 
perturbations set up in a contracting phase can be suc- 
cessfully mapped into a post Big-Bang expanding phase 
(for a more detailed discussion in the ekpyrotic setting 
see e.g. [H, [Mill). 



IX. CONCLUSIONS 

We considered general non-adiabatic cosmological per- 
turbations, allowing for varying w and c 2 es . In this set- 
ting we derived the second-order action and its curvature 
variables in the same generality as has been achieved 
for purely adiabatic modes before. This showed how 
non-adiabatic perturbations give rise to a new conserved 
super-horizon quantity 9. Relating it to the adiabatically 
conserved £ aa - draws attention to potential consequences 
such as enlarged non-Gaussianities. 

Applying the new formalism to an extension of the 
specific background model discussed in [29], we then 
produced both an adiabatic and a non-adiabatic scale- 
invariant solution for structure formation. The specific 
model considered has to be implemented in a contracting 
phase in order to resolve the horizon problem. In this 
context we considered the issue of growing and decaying 
modes. We show that, untypical for contracting models, 
the C\ mode dominant in the corresponding expanding 
case also stays dominant here. 

It will be an interesting task for the future to investi- 
gate whether models with other equations of state and 
speed of sound profiles can produce further expanding 
and (near) scale invariant scenarios for structure forma- 
tion. We hope the formalism developed will be useful in 
this enterprise. 
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X. APPENDIX A: CONSTRUCTING A 
SCHEMATIC ACTION 

Consider a system with equations of motion: 
( + A( V )t = B( V )( C' + C( V )( = £>(»?)€, (74) 

where 77 is a "time" associated with the system, ' denotes 
differentiation w.r.t. 77 and (£,C) are the variables of the 
system. We want to find an action which, when varied 
w.r.t. £ and £, reproduces these equations of motion. 
For this consider the following toy Lagrangian 

(75) 

where O = 0(A) is some time- independent operator. 
The Euler-Lagrange equations then allow us to fix the 
a(rj)..a(r]) terms. Furthermore one can eliminate the 
A{rj) and C(rj) terms by employing substitution schemes 
in analogy with the way we defined the /-measure (|24[) . 
Without loss of generality and with an eye on equations 
((25)) and |26|) we therefore set A(rj) and C(rj) to zero here. 
In fact this procedure gives us a simple tool to combine 
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e.g. terms linear in a variable £ with terms linear in £ 
via a change of variables. After some algebra we then 
obtain the following Lagrangian 

L = ^D(rj)t6t - tyB( v )(d( - K^OC , (76) 

where K± is a constant of integration and we have elimi- 
nated an extra degree of freedom by choosing a — —K\. 
Replacing £ via (|26[) we can therefore finally write down 
the action for the full system: 

S(C,C',V) = J z 2 ((' 2 + B( V )D(r ] K 2 )dr 1 d 3 x, (77) 

where we have defined 

- 6 (78) 



D{ri) 

We are now in a position to introduce the variable 
v = zQ, for which the action (up to a total derivative) 
takes on the following form 

S = J {v 2 + B(r))D(r])v 2 + ^-v^drjifx. (79) 

Varying this action w.r.t. v we obtain the equation of 
motion: 



v - BOrADMv v = 0. (80) 

z 



XI. APPENDIX B : REGIMES OF THE 
EXAMPLE MODEL 

A. A constant equation of state 

Constant equations of state abound in the literature. 
They are invoked e. g. i n the standard de Sitter solution 
or by cyclic models |38| where a constant, high w is pos- 
tulated. In our setup they correspond to /5 = 0. We 
therefore get 



1 + w = 7. 



(81) 



where 7 is a constant. If 7 = 0, the familiar exponentially 
inflating de Sitter solution is reproduced 



a(t) cx e 



Kt 



(82) 



If 7 > 0, then the expression for the scale factor a be- 
comes 



i(t) = aot 3 ~t . 



(83) 



Consequently we are in an inflationary regime ii — 
1) > 0. In particular we have power law inflation for 
7 < 2/3. This model therefore either always shows ac- 
celerated expansion or it never does at all, depending on 
what the constant value of w is tuned to. 



The energy density is now given by 



3 7 2 



(84) 



i.e. we obtain the well-known "scaling solution" for e. At 
t = the scale factor goes to zero and the energy density 
diverges, i.e. we have a Big Bang singularity. 

If a is positive and non-zero, we obtain the following 
behavior of fc-modes 



(k 2 Y oc (i 



-4a-2+- 



< 0. 



(85) 



We therefore have both expanding and contracting so- 
lutions depending on how parameters are chosen. If 
4a + 2 > j- we are in a contracting scenario. Conversely, 
if the inequality is not satisfied, an expanding model is 
the result. However, in both cases scale invariance fails 
to obtain by the same arguments as given in the a > /3 
case. For we get an exponential contribution via the / 
dependence of z 2 IpIS)) . 

In the constant c 2 s limit, on the other hand, (i.e. a = 
0) the condition for an expanding solution interestingly 
becomes 7 < 2/3, so that we have expansion if and only 
if the model displays power law inflation. Such models 
are well known, so we will not discuss them further here. 

It is worth pointing out that the considerations for 
(3 = also apply to generic j3 < models. Remember we 
are considering an equation of state of the form 



W + Wq = 76 



2d 



(86) 



Combining this with the Friedmann equations we get 



+ V3ei(7e 2 / 3 -( U ; -l)). 



(87) 



If wo 7^ 1, the e 2/3 term is suppressed for (3 < 0. Hence 
we again have a scaling solution for the energy density 
e oc t~ 2 . Having gotten rid of a constant of integration by 
shifting the time coordinate, this takes on the following 
form 



e = -(wo - iy 2 r 



(88) 



For a — we are in the power law inflation regime, other- 
wise this results in an exponentially divergent spectrum. 
This case is therefore completely analogous to that of a 
constant equation of state considered above. 



B. Tuned f3 < scenarios 

The cases where /3 < and the equation of state is 
tuned to reproduce the standard perfect fluid equation 
of state, i.e. choosing Wq = 1, are noticeably different. 
The corresponding adiabatic background solution is dis- 
cussed at length in [3l|. Essentially we reproduce the 
expressions presented in IVII Al unless /3 = — j. As a 
general feature such models are naturally drawn towards 
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W = — 1 either as t — >• or for t >> 1. The evolution is 
therefore either set up in an inflationary de Sitter phase 
or eventually exits into one. 

If > ft > — i , the constant a assumption breaks down 
at first sight. However, recall that 

a = a Q e K ^. (89) 

K a is now negative and a is consequently strongly sup- 
pressed, i.e. a << a>o at early times. With respect to the 
background solution this changes very little though. We 
recover / // oc 3Hc s 2 in analogy with the a > ft case. 
Consequently no scale-invariant solution is possible just 
as in the a > ft case, for z 2 is exponentially divergent. 
If ft < —j the background behaves analogously to the 



ft > a > case, a is therefore loitering as a function of 
time. Consequently we can replicate the expression for 
horizon-crossing modes 

(i 2 )'«(r 2+ w)'<o. (90) 

Moving to negative ft has not affected the fact that the 
solution is contracting. On the contrary the contraction 
is more rapid now. The constant c 2 es solution can still be 
implemented and we have identified the negative ft ana- 
logue of the non-adiabatic scale-invariant solution found 
in the positive ft domain. Importantly the consideration 
of growing and decaying modes for the ft > a > case in 
IVIIII also carries over and we find a dominant C\ mode 
again. 
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